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, ABSTRACT 

' In this paper, we investigate the star formation and chemical evolution of damped 

Lyman-a systems (DLAs) based on the disc galaxy formation model which is devel- 
oped by Mo, Mao & White. We propose that the DLAs are the central galaxies of less 
J> ' massive dark haloes present at redshifts z ^ and they should inhabit haloes of mod- 

lO i erately low circular velocity. The empirical Schmidt law of star formation rates, and 

OO . closed box model of chemical evolution that an approximation known as instantaneous 

\l ' recycling is assumed, are adopted. In our models, when the predicted distribution of 

, metallicity for DLAs is calculated, two cases are considered. One is that, using the 

closed box model, empirical Schmidt law and star formation time, the distribution of 
metallicity can be directly calculated. The other is that, when the simple gravitational 
instability of a thin isothermal gas disc as first discussed by Toomre is considered, the 
i-S^ ' star formation occurs only in the region where the surface density of gas satisfies the 

critical value, not everywhere of a gas disc. In this case, we first obtain the region 
Q . where the star formation can occur by assuming that the disc has a flat rotation curve 

$-H ' and rotational velocity is equal to the circular velocity of the surrounding dark matter 

\ halo, and then calculate the metallicity distribution as case one. We assume that star 

d • formation in each DLA lasts for a period of 1 Gyr from redshifts z — 2>. There is only 

one output parameter in our models, i.e. the stellar yield, which relates to the time 
of star formation history and is obtained by normalizing the predicted distribution of 
metallicity to the mean value of 1/13 Zq as presented by Pettini et al. The predicted 
Jh \ metallicity distribution is consistent with the current (rather limited) observational 

data. A random distribution of galactic discs is taken into account. 

Key words: galaxies: formation - galaxies: evolution - galaxies: stellar content. 



1 INTRODUCTION 

Damped Lyman-a systems (DLAs) are high column density 
(the HI column density is higher than 10^'' '^/cm^) absorbers 
detected in the optical spectra of quasars up to relatively 
high redshift (up to z ~ 5), and can be used to probe the 
cosmic chemical evolution, i.e. the globle evolution of gas, 
metallicity and star formation rates in the Universe. DLAs 
are now generally believed to be the progenitors of present- 
day galactic discs (Wolfe et al. 1986; Wolfe 1988). This view 
is supported by the arguments based on mass estimates, 
which are from the column densities and absorber sizes at 
redshifts z ~ 2 — 3, as well as by the kinematic features 
consistent with models of fast-rotating (at velocities of or- 
der 200 kms"^), thick discs (Prochaska & Wolfe 1997a,b, 



1998; Wolfe & Prochaska 1998). At the same time, on the 
basis of [a/Fe] and [N/0] abundance ratios and their large 
scatter, an origin of DLAs in dwarf and low surface bright- 
ness galaxies is also presented (Matteucci, Molaro & Vladilo 
1997; Vladilo 1998; Jimenez, Bowen & Matteucci 1999). Be- 
cause DLAs contain most of the cool, neutral gas in the 
universe at high redshift, they play an important role in 
studies of the galaxy formation, gas process, star formation 
and metal enrichment as well as the large scale structure. 
The pioneering work in this field was done by Meyer, Welty 
& York (1989) and Pettini, Boksenberg & Hunstead (1990). 
In the past few years, Pettini et al. (1994, 1997a) have de- 
voted great efforts in observation of DLAs' column densities 
and metal abundances. They found that, with the increase 
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of observed data, DLAs are generally metal poor with a 
mean value only of 1/13 Zq and are very young galaxies 
in the early stage of chemical enrichment, and the metallic- 
ity distribution of DLAs shows a large spread of nearly two 
magnitudes. They also pointed out that the abundance of 
Cr is relatively poorer than that of Zn due to the depletion 
of dust. 

Lu et al. (1996) used the powerful Keck 10 m telescope 
to study the metal abundances of many heavy elements, 
and presented a very detailed analysis of DLAs including 
comparing their distribution of metallicities, electron den- 
sity and kinematics with those of our Milky Way galaxy, 
etc. Boisse et al. (1998) presented HST spectra for a sample 
of six DLAs with intermediate redshift [za < 1), and con- 
cluded that the available observations may be biased against 
dust-rich absorbers and dust extinction causes a preferential 
selection of QSOs with intervening gas relatively poor in 
metals, dust and molecules. They also strongly suspect that 
surveys of DLAs down to fainter QSOs would reveal more 
systems with large HI column densities and high metallici- 
ties. M0ller & Warren (1998) used HST images to estimate 
the cross-section-weighted mean radius of DLA absorbers at 
high redshift {z > 2), that is they are smaller than 23.6 
kpc and 12.7 kpc for go ~ 0.0 and 0.5, respectively. 
They also found that for qo = 0.5 the space density of DLAs 
at high redshift is more than five times that of the spiral 
galaxies locally, but the problem will be ruled out if qo de- 
creases. The observations of Ha emission done by Bunker 
et al. (1999) and near-77? and IR done by Bechtold et al. 
(1998) also supported the conclusions of M0ller & Warren 
(1998) above. 

Untill now, most of the researches in DLAs are for ob- 
servational constrains. There has, however, been relatively 
little work on its simple prescription. It is worthy of em- 
phasizing that the observation of HST (M0ller & Warren 
1998) strongly suggests DLAs and LBGs may commonly 
associated with each other. Based on their successful disc 
galaxy formation model. Mo, Mao & White (1999) have in- 
vestigated the global properties of LBGs, such as the sizes, 
luminosities, kinematics and star formation rates. They pre- 
sented that the number of density and clustering properties 
of LBGs are consistent with them being the central galaxies 
of the most massive dark haloes present at 2 ~ 3. 

In this paper, under an assumption that the DLAs are 
the central galaxies of less massive dark haloes present at 
2 ~ 3, and they should inhabit haloes of moderately low 
circular velocity, we study how star formation and chemi- 
cal enrichment may have proceeded in them by using the 
simplest models possible with the smallest number of pa- 
rameters. 

The outlines of this paper are as follows. In Section 2, we 
give a simple description of our models; The predicted distri- 
butions of metallicity and neutral hydrogen column density 
for DLAs are presented in Section 3, and some correlations 
will also be shown in this section; The summary and discus- 
sion are presented in Section 4. 



2 MODELLING DAMPED LYMAN-a 
SYSTEMS 

We model the assembly of DLAs in the context of the stan- 
dard hierarchical picture (White & Rees 1978; White & 
Frenk 1991). Structure growth in this model is specified 
by the parameters of the background cosmology and by 
the power spectrum of initial density fluctuations. As an 
illustration, we show theoretical results for a CDM model 
with cosmological density parameter fio = 0.3, cosmolog- 
ical constant f^A — 0.7. The power spectrum is assumed 
to be that given in Bardeen et al. (1986), with shape pa- 
rameter r — 0.2 and with normalization erg = 1.0. We de- 
note the mass fraction in baryons by /b — Qb/^^o, where 
JIb is the cosmic baryonic density parameter. According to 
the cosmic nucleosynthesis, the currently favoured value of 
Qb is SIb ~ 0.019^"^ (Buries & Tytler 1998), where h is 
the present Hubble constant in units of 100 kms~^Mpc~^, 
and so /b ~ 0.063/i~^. Whenever a numerical value of h is 
needed, we take h — 0.7. 

2.1 Galaxy formation 

There are some traditional methods to investigate the galaxy 
formation and evolution. One can, for instance, assume 
the redshift of galaxy formation and star formation his- 
tory, based on the stellar IMF and cosmological model, to 
trace the observations (Tinsley 1980; Pozzetti, Bruzual & 
Zamorani 1996; Jimenez, Padoan & Matteucci 1998). White 
and his collaborators (White & Frenk 1991; Kauffmann & 
White 1993; Kauffmann, White & Guiderdoni 1993; Kauff- 
mann et al. 1999) construct hierarchical models of galaxy 
formation and evolution, and are successful in many obser- 
vations such as the statistical properties of local and distant 
galaxies, etc. Madau, Pozzetti & Dickinson (1998) developed 
a method which focuses on the integrated light of galaxy 
populations as a whole to match the properties of the star 
formation history for the field galaxies. 

In the present paper, we use the Mo, Mao & White 
(1998) (hereafter, MMW) galaxy formation model in assum- 
ing that central galaxies form when collapse of the proto- 
galactic gas is halted either by its angular momentum or by 
fragmentation as it becomes self-gravitating. This model is 
in hierarchical cosmogonies, in which a disc is a thin and 
centrifugally supported structure with an exponential sur- 
face density profile. LBGs can be well understood in MMW 
with large circular velocity Vc and small spin parameter A 
which correspond to compact objects with massive haloes 
(Mo et al. 1999). In this paper, we assume that DLAs are 
the central galaxies of less massive dark haloes present at 
redshifts z ~ 3, and they should inhabit haloes of moder- 
ately low circular velocity. 

According to MMW, if the mass proflle of the disc is 
taken to be exponential 

E(r) = Eoexp(-r/iid), (1) 

then So and Rd that are the central surface density and 
the scale length of the exponential disc, respectively, and 
the rotation curve of the galaxy, are determined uniquely. 
Speciflcally, 

7^..8.8.-kpc(4)(^^^)[^]- (2) 



© RAS, MNRAS 000, 0-?? 



Star formation and chemical evolution of damped Lyman a systems 3 



and 



Eo ~ 380/1 M© pc" 



^0.05 '^^0.05^ 



T4 



250 kms- 



(3) 



where rud is the fraction of disk to halo mass, 14 the circular 
velocity of the halo, A the dimensionless spin parameter, 
H{z) the Hubble constant at redshift z, Hq the present- 
day Hubble constant (See detailed in MMW). Because the 
Hubble constant H{z) increases with redshift, it is expected 
from equations (2) and (3) that the galaxy disc should be 
less massive and smaller but have a higher surface density 
at higher redshift. One can find if K, A and cosmogony are 
given, _Rd, which dominates the disk profile, will not change. 

For a given cosmogony, the mass function of dark mat- 
ter haloes at redshift z can be estimated from the Press- 
Schechter formalism (Press & Schechter 1974): 



diVh(A4,z) = 



2 po 4(z) dlnA(Mh) 
7rA/hA(A/h) dlnAfh 



X exp[ 



(4) 



where 5c{z) — Sc{0){l + z)g{0) / g{z) is the linear overdensity 
corresponding to collapse at redshift z, and g{z) is the linear 
growth factor relative to an Einstein-de Sitter universe at z, 
(5c (0) ~ 1.686, A(Afh) is the rms linear mass fluctuation at 
2; = in a sphere which on the average contains mass Alh 
by A/h ~ {^■7v/3)poR^ , po is the mean mass density of the 
universe at 2: = 0. The relation between the mass Afh and 
circular velocity 14 of a halo is in equation (2) of MMW, 
and a more detailed description of the PS formalism and 
the related issues can be found in the Appendix of MMW. 

Based on the PS formalism and the A distribution, with 
the latter being a log-normal function with mean In A = 
In 0.05 and dispersion (Tin a = 0.5 (MMW), we can generate 
Monte Carlo samples of the halo distribution in the A^h-A 
plane at redshifts z = 3. We can then use the galaxy forma- 
tion model of Section 2.1 to transform the halo population 
into a DLA population. Finally, suppose that each halo at 
z ~ 3 has a central galaxy with a star formation rate (SFR) 
that is a monotonic function of halo mass, and only a negli- 
gible fraction of haloes hosts a second galaxy bright enough 
to be seen. The observed DLAs then correspond to haloes 
at a ~ 3. Here a lower-limit for Vc of 50 km/s of haloes is 
chosen, and a random distribution of galactic discs is taken 
into account. 

It should be pointed out that the interaction between 
discs and bulges is not considered in the present study. This 
interaction will become more and more important for very 
compact objects which correspond to A < 0.025 (MMW). 
Here, we treat galaxies with A < 0.025 as A = 0.025. It can 
be easily estimated from Warren et al. (1992), Cole & Lacey 
(1996) and Steinmetz & Bartelmann (1995) that, there are 
only 10 percent galaxies with A < 0.025. This will not influ- 
ence our main results. 



2.2 Star formation rates 

For star formation rates (SFR), we take the empirical 
Schmidt law (Kennicutt 1998) 



SsFR = a{ 



where 



Sg(r) 



IM0PC- 



'Meyr- 



'1 -2 
pc , 



a = 2.5 X 10" 



,6= 1.4, 



(5) 



(6) 



respectively. Here Esfr is the SFR per unit area, and Sg(r) 
the gas surface density. This star formation law was derived 
by averaging the cool gas density over large areas on spiral 
discs and over starburst regions. 



2.3 Closed box model 

For a closed box model, the metallicity Z of gas is (Binney 
& Merrifield 1997; Tinsley 1980) 



Z-Z,^ -2/ln(Eg(r)/Ego(r)) =yln^- 



(7) 



where Zi is the initial metallicity of cool gas that is taken 
to be 0.01^0 in the present paper, y the stellar yield that is 
obtained by normalizing the predicted distribution of metal- 
licity to the mean value of 1/13 Zq as presented by Pettini 
et al. (1997a) in this study (see next section), Ego(r) the 
initial gas surface density and jj, the gas fraction, and 



dEg(r) 
dt 



(1 — -Rr)EsFR, 



(8) 



where Ri is the returned fraction of stellar mass into the 
ISM, and we take Rr — 0.3 for a Salpeter IMF (see Madau et 
al. 1998). We assume that at initial time {t^ = 0), Esfr — 
and Ego(r) = Eo(r) exp(— r/_Rd), with the latter being given 
by equation (1). After a period of T, the surface density of 
the survived cool gas can be obtained by equations (5) and 
(8), that is 



Eg(r) = Ego(r)(l + r/r)- 



and the current SFR is 



Esfr = aEgo(r)(l +T/t) 



where 



a{b- 1)(1 - R, 



(9) 



(10) 



(11) 



As we know, in the simple closed-box model of chemical 
evolution, all metals produced by stars are retained within 
the galaxy. The metallicity of the galaxy is simply deter- 
mined by the stellar yield y, deflned as the mass of metals 
produced per solar mass of long-lived stars that are formed 
(Tinsley 1980). It is worthy to note that, the instantaneous 
recycling approximation is not valid as regards the contri- 
butions to metal enrichment of either type la supernovae or 
of stars of a few solar masses (Binney & Merrifield, 1997). 
Pagel (1997) has refined to take into account the time de- 
lay between star formation and the ejection into the ISM of 
heavy elements. 

In this paper, we use [X/H]^ for the logarithm metal- 
licity, and all the histograms are normalized to unity. 



* We use the normal notation where [X/V] = log(X/Y) 
log(X/Y)Q. 
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2.4 Criterion for star formation 



According to Toomre (1964), in the case of a gaseous disc, 
instability for star formation is expected if the surface den- 
sity of gas exceeds the following critical value 



3.36G' 



(12) 



where a is taken to be 2 (van der Kruit & de Grijs, 1999); 
and 



V r AV 
1.41^(1 + :^^) ; 
r V fir 



(13) 



and (Tg is a gas velocity dispersion, we take it to be 6kms~^ 
(Kaulfmann 1996), which is consistent with the observa- 
tional results at any redshift till now (see Larson 1998). 

When this criterion is considered, the star formation 
occurs only in the region where the surface density of gas 
satisfies the critical value, not everywhere of the gas disc. 



3 RESULTS 

In this section, we present the predicted normalized distri- 
butions of metallicity and HI column density for the DLA 
population. As explained by Pettini, King & Smith (1997b), 
[Zn/H] is a straight forward measure of the degree of metal 
enrichment analogous to the stellar [Fe/H]. So, we use the 
mean value of [Zn/H] as presented by Pettini et al. (1997a) 
to normalize the predicted distribution of metallicity for 
DLAs in order to obtain the steller yield. A random dis- 
tribution of discs is taken into account. 



3.1 Distribution of metallicity without 
considering the Toomre's criterion 

At first, we do not consider the criterion for star forma- 
tion which is developed by Toomre (1964). Based on the 
Monte Carlo samples at z = 3, we can calculate the metal- 
licity Z of DLA population by the formulas from (5) to (11) 
by assuming that, at initial time {ti, — 0), Ssfr — and 
^go(^) ~ So exp(— r/i?d); and the star formation time T 
is 1 Gyr. Eo and Rd are decided by equations (2) and (3). 
The stellar yield is obtained by normalizing the predicted 
distribution of metallicity to the mean value of 1/13 Zq 
as presented by Pettini (1997a). It is 0.22 Zq. The results 
are plotted in Fig. 1. For comparison, the metallicities of 
18 DLAs at redshifts z ~ 3 (2.3 < z < 3.4) that were ob- 
tained by Meyer, Welty & York (1989) and Pettini et al. 
(1994, 1997a) |^ are also presented. This figure shows that, 
the observational data that have a spread of nearly two or- 
ders of magnitude agree with our model, but the predicted 
metallicity distribution cannot spread as widely as the ob- 
servations. The observational data are still quite sparse, and 
larger and more complete samples are needed. In order to 
test whether our model fits the abundances of DLAs, we 
provide the Kolmogorov-Smirnov Test. The results are that 
the maximum value of the absolute difference is 0.33, and 
the significance level probability is 69.94 percent. 

t The solar system abundance of Zn is from the compilation by 
Anders & Grevesse (1989), i.e. log(Zn/H)0 = --7.35. 
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Figure 1. The distribution of metallicities for DLAs. The solid 
histogram gives the model prediction and the dashed histogram 
shows the observational data in Meyer et al. (1989) and Pettini 
et al. (1994, 1997a). 



3.2 Distribution of metallicity with considering 
the Toomre's criterion 

When we take into account the Toomre instability criterion 
of star formation, the star formation occurs only in the re- 
gion where the surface density of gas must be larger than 
Ec(r-). So, we first obtain the regions where the star forma- 
tion can occur based on the formulas (1), (12) and (13), as- 
suming that the disc has a fiat rotation curve and rotational 
velocity is equal to the circular velocity of the surrounding 
dark matter halo. If the observation is not in these regions, 
the predicted distribution of metallicity for DLA population 
is initial one of the cool gas, i.e. 0.01 Zq. When the obser- 
vation is just in these regions, the metallicity of the DLA 
population can be obtained as section 3.1. In this case, the 
stellar yield is 0.70 Zq. Fig. 2 presents the normalized model 
distribution of metallicity. For comparison, the metallicities 
of 18 DLAs at redshifts z ~ 3 (2.3 < z < 3.4) that were 
obtained by Meyer, Welty & York (1989) and Pettini et al. 
(1994, 1997a) are also plotted. This figure shows that, the 
model data which spreads more than two magnitudes can 
match the observed ones (Meyer et al. 1989, Pettini et al. 
1997a), although the model is a little metal rich. As in Fig. 1, 
we also provide the Kolmogorov-Smirnov Test. The results 
are that the maximum value of the absolute difi'erence is 
0.25, and the significance level probability is 84.75 percent. 



3.3 Distribution of HI column density 

HI column density is the major characteristics of DLAs, be- 
cause it is the first physical parameter that we can obtain 
from observations of QSO absorption lines. If we assume that 
the survived gas from the star formation in a disc is dom- 
inated by the neutral hydrogen, the model column density 
of neutral hydrogen can be obtained by the formulas from 
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Figure 2. The distribution of metallicities for DLAs. Tlie solid 
histogram gives the model prediction and the dashed histogram 
shows the observational data in Meyer et al. (1989) and Pettini 
et al. (1994, 1997a). 



Figure 3. The distribution of HI column densities for DLAs. 
The solid histogram gives the model prediction and the dashed 
histogram shows the observational data in Meyer et al. (1989) 
and Pettini et al. (1994, 1997a). 



(5) to (11). The predicted distribution of HI column den- 
sity for DLAs without and with considering the Toomre's 
criterion are shown in Figs. 3 and 4, respectively. For com- 
parison, the HI column densities of 18 DLAs at redshifts 
z ^2, (2.3 <z < 3.4) that were obtained by Meyer, Wehy & CD 
York (1989) and Pettini et al. (1994, 1997a) are also plotted. ^ 
These two figures present that, the observational neutral hy- ^ 
drogen column densities spread narrower than the predicted r-j 
ones. The reason is, perhaps, that the extinction may results ^ 
in a preferential selection of QSOs with systems displaying 
a low HI and H2 content. 1) 

In fact, Boisse (1995) has already noticed that, the ab- ^ 
sence of very large N(HI) values (e.g. > 10^^ cm~^) which +-> 
would be quite easily recognized in the numerous low reso- ^ 
lution optical spectra obtained so far, supports the reality qj 
of cut-off caused by extinction. 

Surveys of DLAs samples drawn from the observation 
of fainter QSOs may reveal more systems with large N(HI), 
and spectroscopic surveys of fainter QSOs can lead to a more 
representative view of the whole diversity of the DLAs. 




20 



21 22 
Log(N(HI)) 
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3.4 Relation between metallicity and HI column 
density 

In order to present the predicted relation between metallicity 
and HI column density for DLAs, we randomly select 100 
samples from the DLA population by their cross-sections 
because the observations of DLAs are based on the QSO 
absorption lines. 

In Figs. 5 and 6, we plot the predicted distribution of 
[Zn/H] as a function of HI column density for DLAs without 
and with considering the Toomre's criterion, respectively. 
For comparison, the observational data of 18 DLAs at red- 
shifts 2; ~ 3 (2.3 < z < 3.4) that were obtained by Meyer, 
Welty & York (1989) and Pettini et al. (1994, 1997a) are 



Figure 4. The distribution of HI column densities for DLAs. 
The solid histogram gives the model prediction and the dashed 
histogram shows the observational data in Meyer et al. (1989) 
and Pettini et al. (1994, 1997a). 



also plotted. The predicted correlations (especially, without 
considering the Toomre's criterion) are that, the higher the 
column density, richer the metallicity. In fact, Boisse et al. 
(1998) have already emphasized that, in an unbiased sample, 
the large N(HI) values should correspond to the innermost 
parts of galaxies where the metallicity is presumable higher. 
The observational data show a different trend, i.e. the metal- 
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Figure 5. The relation between metallicity and HI column den- 
sity for DLA population. The blank circles and crosses represent 
the model prediction and the observational data in Meyer et al. 
(1989) and Pettini et al. (1994, 1997a), respectively. 
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Figure 6. The relation between the metallicity and HI column 
density, where blank circles and crosses represent the model pre- 
diction the observational data in Meyer, Welty & York (1989) and 
Pettini et al. (1994, 1997a), respectively. 



White (1998). One is without considering the Toomre's cri- 
terion, another is with considering the Toomre's criterion. 
The star formation rate is chosen to be in the form of the 
empirical Schmidt law (Kennicutt 1998). The approximation 
of instantaneous recycling is assumed. The predicted distri- 
bution of metallicity for DLAs spreads larger than one order 
of magnitude, which is consistent with the observations by 
Meyer, Welty & York (1989) and Pettini et al. (1994, 1997a). 

Although the star formation time of 1 Gyr, which we 
adopted, is arbitrary, its order is consistent with that of 
observations (Bechtold et al. 1998). It is no problem for us 
to choose it as a typical time of star formation to investigate 
the global properties of DLAs. 

In our models, the star formation time scale in the disc 
is implicitly assumed to be the same everywhere. This will 
lead to the slow consumption of gas in the central parts of 
the disc. In fact, as we know well, the star formation in the 
central region is very rapid, and the central region should 
correspond to higher metallicities and lower HI column den- 
sities. Furthermore, we treat all the survived gases as HI. 
There actually exist H2, ionized gases and other elements. 
This is also leads to a overestimation of HI column density, 
especially in the central region of the disc. In addition, the 
DLAs that correspond to have small Vc, should form earlier 
and the merging could take place more frequently. So, the 
DLAs are to a certain extent more complicated. 

There is only one output parameter in our models, i.e. 
the stellar yield, which relates to the star formation time 
scale and is obtained from the normalization of the predicted 
distribution of metallicity to the mean value of 1/13 Zq as 
presented by Pettini et al. (1997a). The metallicity of the 
initial cool gas is assumed to be of 0.01 Zq. According to 
Pagel (1987), the stellar yield is 0.3 Zq for the disc clusters 
and 0.4 Zq for the solar neighborhood. The stellar yields in 
our models are also quite similar to that of LMC of some 
0.2 Zq (Binney & Tremaine 1987). 

Finally, we must point out, although our models are 
very simple without many physical processes being consid- 
ered, such as the interaction between discs and bulges, merg- 
ing effects, gas reheating and supernovae feedback, etc., it 
can predict some properties of DLAs. 
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licity is richer with the decreasing of column density because 
of the extinction bias as pointed out by Boisse et al. (1998). 



4 SUMMARY AND DISCUSSION 

In the present paper, we establish two simple scenarios to 
investigate the nature of DLAs at redshifts z ~ 3 based on 
the disc galaxy formation model developed by Mo, Mao & 
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